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On the basis of a previous tiieoretical approach to the plastic flow of highly refined materials, a 
physical explanation for diffusion bonding is essayed, which yields closed-form equations relating 
the bonding progress with time, temperature, applied pressure and the constants characterizing the 
material. Excellent agreement with experiment is attained, with no adjustable parameter. In the 
novel scheme, diffusion bonding is caused by the interpenetration of the two sufaces at the grain 
level. The process is driven by the strong tensile stress field induced in the plane of the interface by 
the plastic deformation in the normal direction. The grain boundaries of each joining surface yield 
to host grains of the other surface, releasing this way the internally generated tensile stresses. Voids 
gradually close with the increment of the interpenetrated areas. In this scheme bonding is not a 
matter of contacting and atomic interdiffusion, but of grain exchange. 

PACS numbers: 62.20.fq, 62.20.fk, 62.20.mq, 83.50.Uv 



Diffusion bonding (DB) is a solid state high tempera- 
ture process for making a monolithic joint through mat- 
ter transportation across the interface between the ma- 
terials being joined, with no melting at the bond line. It 
is a costly time-consuming technique for producing top 
quality bonds, so that neither metallurgical discontinu- 
ities nor voids can be detected all over the former in- 
terface. Though not worthy for mass production assem- 
bly lines yet, millions people are benefited every day by 
diffusion bonded artifacts. The titanium hollow blades, 
their internal reinforcing honeycomb structure, and the 
central disc of the wide chord fan of a present day turbo- 
fan commercial aeroengine, which generates about 80% 
of the total engine thrust at take-off, integrate a sin- 
gle component diffusion bonded unit. On the basis of 
a recently published theoretical approach to the plastic 
flow of highly refined materials, we essay here a physi- 
cal explanation for DB and derive a closed-form equa- 
tion relating bonding time, temperature, applied pres- 
sure and the constants characterising the material, which 
exhibits excellent agreement with available experimental 
data, and with no adjustable parameter. In our theo- 
retical approach, DB has little connection, if any, with 
atomic diffusion. 

For two decades, DB combined with superplastic form- 
ing (SPF) has been a standard technique for producing 
aircraft structural components and engine parts, partic- 
ularly because SPF yields practically no residual stress 
and neither spring back from die after releasing the form- 
ing pressure. Besides, DB eliminates fasteners, yielding 
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stronger more reliable unions, saving weigth and allow- 
ing better design [2 [5]. However, the suitability of the 
process for mass production is yet questionable, mainly 
by the long bonding times involved and because DB does 
work just for a short list of known materials. The precise 
knowledge of the microscopic mechanisms taking part in 
the process, leading to the physical equations governing 
its evolution, would be of great value as scientific knowl- 
edge, as well as a tool for the search of new improved 
materials and process optimisation. 

The theoretical literature on the subject assumes that 
any pair of clean metallic surfaces will bond if they are 
brought together within the range of interatomic forces. 
Bonding happens by atomic interdiffusion across the in- 
terfaces at the contacting surface sectors, which initially 
are isolated and scarce by the surface roughness. Con- 
tacting area gradually increases by plastic deformation 
and atomic interdiffusion. The large differences in bond- 
ing ability observed for the various alloys are ascribed to 
the chemical and mechanical stability of their surface ox- 
ide layers, which may obstruct atomic exchange [HISlllj. 

In our scheme the process is quite different. It takes 
place mainly by the active migration of grains that cross 
the joint line in the two senses, rather than intergrain 
diffusion of individual atoms. Grain exchange is much 
more efhcient and demands surface proximity of the or- 
der of just the grain size, instead of the much smaller 
atomic distances. The driving force for the interfacial 
grain exchange arises from the internal stress which plas- 
tic deformation necessarily induces in the plane normal 
to the flow direction. 

DB has been observed in highly refined materials con- 
stituted by small equiaxed grains, whose plastic deforma- 
tion conserves the grain mean size d and overall shapes. 



Thus, grains slide past each other over relatively long dis- 
tances, accommodating their shapes to preserve matter 
continuity, with grain boundaries maintaining mechani- 
cal integrity and coherence. Hence grains evolve by the 
action of two kinds of forces that may be of very different 
strength: those making their boundaries to slide and the 
induced local forces responsible for the continuous grain 
reshaping. The stronger one will determine the plastic 
flow rate. 

Recent papers show that excellent agreement with ex- 
periment is attained from ascribing the rate control to the 
intergrain forces that make adjacent grains to slide [SHS]. 
These forces have been examined by theory [3 [6] , Mon- 
tecarlo computer simulation |11| . experiment |10| and 
molecular dynamics simulation \p. Said studies clearly 
show that the onset of grain boundary sliding is associ- 
ated to a boundary mechanical instability which demands 
that the shear stress in the plane of the common bound- 
ary of the sliding grains be higher than a critical value 
Tc- Then the grains start sliding with a relative velocity 
Av which may be proportional to: (a) The difference 
between the in-plane ovcrcritical shear stress and the 
threshold stress Tc- (b) The bare in-plane shear stress, 
and Aw jumps from zero to the proper value when Tc is 
surpassed. We refer to these two options as force models 
A and B, respectively [5]. 

As the two models can be fully worked out without 
additional assumptions, comparison with experiment de- 
termines what model is the right one for a specific sys- 
tem. It has been shown that model A works for a series 
of aluminium and titanium alloys [6l [7] , fitting the exper- 
imental data within the experimental uncertainties, and 
model B gives very good results for the steel Avesta 2304 

i- 

The technical details for deriving the plastic flow of 
the polycrystalline medium from the hypotheses stated 
above are in the literature [SHH] . Summarising, given the 
stress tensor (aij), i, j = x, y, z, the force law A or B de- 
termines the relative velocities Av{9^ (j)) between adjacent 
grains. Angles and (j) characterise the grain boundary 
orientation. One can readily realise that 
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where v is the velocity field of the material medium and 
symbol {■ ■ ■)j stands for the average over all boundary 
orientations along the Xj axis in the positive sense. (No- 
tice that we denote either i,j = x,y,z or Xi = x,y,z, 
i — 1,2,3). Making the averages in explicit way, closed- 
form equations relating dvi/dxj with the stresses are de- 
rived from Eq. (11]). Rotation invariants like V • v and 
V X ■(;, or the components iij of the strain rate tensor, 
which is the symmetric part of tensor (fTl), can be con- 
structed with them. 

For uniaxial external stress a along the z axis and 
cylindrical symmetry with respect to the same axis, so 



that <Txx = o'yy = a± and (Tij — for i ^ j, the strain 
rate along z is shown to be 
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where s is the sign of the deviatoric stress cr -\- p, being p 
the pressure p = —{a + 2crj_)/3, and the coefficient Q{p) 
is the proportionality factor between the relative speed of 
adjacent grains and the resolved ovcrcritical shear stress. 
Temperature T apart, Q(p) depends only on the pres- 
sure invariant p because it must keep unchanged for any 
boundary orientation. The auxiliary variable 9, is de- 
fined by 
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and a takes the values 1 or for force models A or B, 
respectively. 

Eq. ([3]) discloses a remarkable consequence of the exis- 
tence of a finite threshold stress for grain sliding. As 
V • z/ = V/V ^ , where V/V is the dilation rate 
per unit volume, the specific volume is not conserved 
in plastic flow if Tc 7^ 0. As grain volume variations 
can only be elastic, one can recall Hooke's law and write 
W ■ V ~ —{1/B)p, where B is the bulk elastic modulus. 
Combining this with equation (IS]) it follows the equation 



Bt,Q{p) 
' 2d 



1 - cos(26'c 
sin(20,) 



261, 



7rsin(26'c) 



+ -{\-2a)cos{29c) 
■n 

(5) 

for the evolution of p on plastic flow. 

Pressure p monotonically increases or decreases for 

plastic axial stretching (s > 0) or shrinking (s < 0), 

respectively. As a is the externally applied stress, the 

increase of \p\ must be attributed to the variation of 



C^-L = -2(^ + 3^), 
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which represent an internal radial stress field normal to 
the plastic flow direction. Settling (Tj_ = (or p = — cr/3) 
at start, when the strain e = 0, is a natural initial con- 
dition. However, as the deformation proceeds, a±^ be- 
comes finite, taking negative values on stretching, and 



hence contributing to neck formation, and positive ones 
on axial compression. Both p and a± grow at a high rate 
because B ^ 10^ MPa. If the plastic flow is interrupted, 
the deformation induced stress field (t± will remain as a 
residual stress. 

We must remark that the equations written above as- 
sume cylindrical symmetry. In strict rigor, they hold for 
an axially symmetric sample or at any point of an infinite 
medium. Symmetry breaking surfaces make the problem 
much more complex. However, there is a practical sit- 
uation for which Eqs. (pi) and (Is]) hold strictly: if the 
material is finite in the z direction, but infinite in the xy 
plane, one can anyway assume cylindrical symmetry at 
any point of it. Hence, Eqs. (pi) and (Is]) are obeyed by an 
extended flat plate subjected to a compressive stress a 
normal to its plane, and strong enough to produce plas- 
tic deformation. In such case the material has no means 
to relax the tensile stress a±, which progressively builds 
up in the xy plane upon plastic compression in the z 
direction. 

However, there is again a practical situation for which 
the extended plate is able to relax the induced tensile 
stress, at least partially, at the surfaces. If we compress 
together two plates with surface irregularities of the or- 
der of the grain size d, the grains of one of them may cede 
to a± and open up if the consequent interstice is simul- 
taneously filled by one or more grains provided by the 
other surface. The process is reciprocal, and the inter- 
penetration of the two surfaces at the grain level provides 
a mechanism to release the tensile cumulative stress a± . 
The interpenetration of the contacting surfaces makes 
their areas to increase, filling this way the voids left by 
the surface irregularities. 

The rate of bonding is dictated by the growing rate of 
the total interpenetrating surface area A. A good cri- 
terion is to assume that A varies at a speed that keeps 
crj_ = 0. We have that A/ A = V • w - dvjdz. Hence A/ A 
is obtained substracting Eq. ([2| from Eq. (|3|, with s — 
— 1 and p — — cr/3, or, equivalently, sin(20^ = 2tc/\(j\. 
Additionally, the coefficient Q{p) was studied in Refs. 5 
and 6 The mechanisms for grain boundary sliding con- 
sidered there give 
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where fcs is the Boltzmann constant, the coefficient Co 
depends only on the grain size d, the constant cq is the 
energy necessary for evaporating a crystal vacancy from 
the grain boundary, and fl* measures the sensitivity of 
this energy to stress. Combining all this, one has that 
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where the subscript of {A/A)q indicates that equation ([8 
gives the initial grow rate of the bonded area fraction, 
when a_\_ — Q. 

The condition crj^ = will cease only once voids be al- 
most filled. Hence, what matters in practice is the initial 
binding rate A, and it can be assumed that A vanishes 
when cTj^ becomes finite because A has no free room to 
grow. On this basis one can write the equation 



A = \^{A,-A)., 
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where Ai is the total area to be bound, Aq is the initial 
contact area, constituted by Ai/Aq small islands which 
are assumed already joined. A is the bonded area at 
time t and Ai — A the unbound area. The coefficient A 
measures the intrinsic ability of the material for making 
DB joints. Eq. (pi) has the solution 
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which gives the time dependent area fraction bonded. 



y 10.0 





steel Avesta 2304 
Force model B 




^^ 


€„ = 2.177 [eV] 




,,00^ 


J==^ 


fi'= 1.077 xlO"^' 


[cm'] _/^^^ 






Co= 2.526 xlO^[ 


''^^ 






^-^ 


Y 


r[K] 


T^, [MPa] 






1173 
1193 
1213 


1.10 
0.85 
0.35 



10 10 ^ 

STRAIN RATE £ [s"''] 



FIG. 1: Flow stress versus strain rate for the steel Avesta 
2304 at r = 1173, 1193 and 1213 K. Circles represent experi- 
mental data taken from Ref. iT7and solid lines depict Eq. H 
with a = and the constants shown in the insets. The pa- 
rameter Tc exhibits some temperature dependence. 

Pilling, Wang and Ridley [12] authored one of the very 
scarce experimental observations of the time progression 
of DB available in the literature. Together with the time 
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FIG. 2: The variation of the area fraction bonded with time 
for steel Avesta 2304 at T = 1213 K and cr = 3.0 MPa. Circles 
represent experimental data taken from Ref. [T^ The solid 
line represents equation ( 10 1 with A given by Eq. ([8| and the 
material constants obtained from the fit shown in Fig. IT] The 
theoretical curve has no adjustable parameter. 



evolution A{t)/Ai of two bonding surfaces of the steel 
Avesta 2304 at T = 1213 K and a = 3.0 MPa, they re- 
ported also the flow stress a for plastic stretching versus 
strain rate e at several temperatures. Therefore, their 
data allows us to test Eq. ([2|, and Eqs. ^ and ([lO| in- 
dependently. Fig. IT] shows the fit given by Eq. (pi) to the 



flow stress data (e « 0) at three temperatures, which 
gives us precise values for the constants eq, Cq, ri* and 
Tf. of the material. 

Inserting these material constants in Eq. Q, together 
with a = 3.0 MPa and T = 1213 K, it is obtained that 
A = 4.11 X IQ-^s-i = 2.47 x lO'^min^^ On the other 
hand. Pilling et al. [T^ reported that the two joining steel 
surfaces in the DB experiment initially display very small 
contact sectors, of less than 4 /um, distant / = 50 — 55 /im 
between them. We take for the ratio Ai/Aq = (l/d)'^ 
because the size of the initial bound spots cannot be 
smaller than the grain size d = 5 /im. Therefore, tak- 
ing / = 52 /im it turns out Ai/Aq — 108.16. The solid 
line in Fig.[2]depicts Eq. ( 10 1 after replacing in it all these 
constants, together with the experimental data of Pilling 
et al. The agreement is quite impressive, particularly 
because it was attained with no adjustable parameter. 

In summary, we obtain close agreement with exper- 
iment from interpreting DB as the interpenetration of 
the two sufaces at the grain level. The process is driven 
by the strong tensile stress field induced in the plane of 
the interface by the plastic deformation in the normal di- 
rection. At each point of contact the grain boundaries of 
one surface yield to host grains of the other surface, and 
reciprocally, releasing this way the internally generated 
tensile stresses. Voids close by the gradual increment of 
the interpenetrated contacting areas. Hence, bonding is 
not a matter of contacting and atomic interdiffusion, but 
of grain exchange. 
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